Determining the density content of symmetry energy and neutron skin: an empirical 
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The density dependence of nuclear symmetry energy remains poorly constrained. Starting from 
precise empirical values of the nuclear volume and surface symmetry energy coefficients and the 
nuclear saturation density, we show how in the ambit of microscopic calculations with different 
energy density functionals, the value of the symmetry energy slope parameter L alongwith that for 
neutron skin can be put in tighter bounds. The value of L is found to be L= 64±5 MeV. For ^°**Pb, 
the neutron skin thickness comes out to be 0.188 ±0.014 fm. Knowing L, the method can be applied 
to predict neutron skins of other nuclei. 
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In recent times, there is a cultivated focus on a better 
understanding of the density properties of the symmetry 
energy of nuclear matter. Particular attention is given to 
constrain in a narrow window the value of the symmetry 
energy slope parameter L at the nuclear matter satura- 
tion density po- In terrestrial context, this parameter af- 
fects the nuclear binding energies [l| and the nuclear drip 
lines and has a crucial role in determining the neutron 
density distribution in neutron-rich nuclei. In astrophysi- 
cal context, it is also of seminal importance. The pressure 
Pn (=3po^) of neutron matter at pq influences the radii 
of cold neutron stars. The cooling of proto-ncutron stars 
through neutrino convection Q , the dynamical evolution 
of the core-collapse of a massive star and the associated 
explosive nucleosynthesis depend sensitively on the sym- 
metry energy slope parameter In the droplet model 
0, Q of the nucleus, the neutron skin is proportional to 
L, a linear correlation between the neutron-skin thickness 
of the nucleus and neutron-star radius 0] could thus be 
envisaged. 

The symmetry energy slope parameter is defined as 



~ dp 



(1) 



where Cy (p) is the volume symmetry energy per nucleon 
of homogeneous nuclear matter at density p. Estimates of 
L are fraught with much uncertainties. Isospin diffusion 
predicts L = 88± 25 MeV [1, |^, nucleon emission ratios 
[lO| favor a value closer to L ~ 55 MeV, isoscaling gives 
L ^ 65 MeV [lljLAnalysis of giant dipole resonance 
(GDR) of 208pb [l3l is suggestive of L - 45-59 MeV, 
whereas pygmy dipole resonance [31 in ^^Ni and ^■^^Sn 
would yield an weighted average in the range L= 64. 8± 
15.7 MeV. Of late, from a sensitive fit of the experimen- 
tal nuclear masses to those obtained in the finite-range 
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droplet model [l[, the value of L could be fixed in the 
bound L = 70 ±15 MeV. Astrophysical observations of 
neutron star masses and radii reportedl y p rovide tighter 
constraints to L to 43 < L < 52 MeV |li| within 68 % 
confidence limits. 

Correlation systematics of nuclear isospin with the 
neutron skin thickness [l^, [l^ for a series of nuclei in the 
framework of the nuclear droplet model has been under- 
taken by the Barcelona Group. This has yielded a value 
of L =75 ±25 MeV. The neutron skins were measured 
from antiprotonic atom experiments [itI . [l8l | , systematic 
uncertainties involving model assumptions to deal with 
strong interaction is therefore unavoidable. The novel 
Pb-radius experiment (PREX) at the Jefferson Labora- 
tory has now been attempted through parity-violation 
in electron scattering as a model-independent probe of 
the neutron density in ^°®Pb [l^. The neutron skin 
Rskin = Rn — Rp was found to be 0.33 ^^'^g, fm, where 
Rn and Rp are the point neutron and proton root-mean 
squared (rms) radii. A reanalysis yielded the value to 
be 0.302 ±0.175 fm The droplet model as weU as 

calculations with class of different interactions, Skyrme 
or relativistic mean- field (RMF), have now clearly estab- 
lished that the neutron skin thickness of ^"^Pb is strongly 
correlated with the density dependence of symmetry en- 
ergy around saturation |2l| - [24j . In the backdrop of this 
information, the large uncertainty in the experimental 
neutron radius of ^"^Pb seems to be of not much help in 
putting L in a tighter bound. 

Nuclear dipole polarizibility cud has been suggested 
[25l.[26j as an alternative observable constraining the neu- 
tron skin. The recent high resolution {p,p') measurement 
[131 of «£) yields the neutron skin thickness of ^°®Pb to be 
0.156^0 021 fn^' I^nt the model dependence [1^ in the cor- 
relation between Rskin and ajj assessed in systematic cal- 
culations in the framework of nuclear density functional 
theory is seen to shift the value of Rskm to 0.168±0.022 
fm. 

In this communication, we suggest a new method for 
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exploiting the empirical information on the volume and 
surface symmetry energy coefficients, Cy{po) and Cs- 
The symmetry coefficient of a finite nucleus asym{-A) can 
be parametrized as 



(2) 



These coefficients have recently been meticulously stud- 
ied [2§| by using the double differences of "experimental" 
symmetry energies. This has the advantage that other 
effects (such as pairing and shell effects) in symmetry 
energy can be well canceled out from the double differ- 
ences for neighbouring nuclei. The correlation between 
the double differences and the mass number of nuclei is 
found to be very compact yielding values of Cy{po) and 
Cs as 32.10 ±0.31 MeV and 58.91 ±1.08 MeV, respec- 
tively. The uncertainties in these symmetry components 
are much smaller than those found earlier. We show be- 
low that these 'experimental' values of C„ and Cs along- 
with empirical information of the proton rms radius in a 
heavy nucleus yields the value of L within narrower lim- 
its; precision information on the neutron skin of nuclei 
also follows from the analysis. 
We start with the ansatz 



a(p) = a(po)(-r, 

Po 



(3) 



where 7 measures the density dependence of the sym- 
metry energy. In a considerable density range around 
Pq this ansatz is found to be very consistent with the 
density dependence obtained from the nuclear equa- 
tion of state (EOS) with different interactions [1, lol. |30| 
and also from experiments in intermediate-energy heavy- 
ion collisions [ifl |il|. At very low densities, however, 
there arc some small deviations. From Eqs. (1) and 
(3), L ~ 3^Cy{pQ) and the symmetry incompressibility 

Ksym = QPo^-§^\po = 3L(7 - 1). One can expand the 
volume symmetry coefficient around po as 



Cy{p) = C^ipa] 
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In the above expansion, keeping terms upto second order 
has been found to be a reliable approximation in our sub- 
sequent calculations. In terms of the symmetry energy 
slope parameter and symmetry incompressibility, Cy{p) 
is then given as. 



C,,{p) = Cy{po)-Le + 



^sym ^2 



(5) 



where e = (po - p)/ (3po)- 

For a finite nucleus of mass number A, the symmetry 
coefficient asym{A) is always less than Cv{po). The co- 
efficient asym{A) can be equated to Ci,(p^) where pA is 
an equivalent density, always less than po- Using rela- 
tions (2) to (5), we show below how pA and 7 (hence L 
and Ksym) can be calculated. From Eqs. (2) and (5) it 
follows that 



r7_ 4-1/3 



if, 



where ca = (po - Pa)/{3po)- From Eq. (6), 

Cs = SCMA'/'bcA - ^7(7 - 1)4]- (7) 

The value of Cs is an empirically determined constant 
(58.91 ±1.08 MeV). 

In the local density approximation, the symmetry co- 
efficient asym{A) can be calculated [30| as 
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'■symy^'-l — j^j^2 



d\p{r)Cy{p{r))[X{r)] 



(8) 



where Xq is the isospin asymmetry {={N — Z)/A) of 
the nucleus, p{r) is the sum of the neutron and proton 
densities inside the nucleus and X{r) is the local isospin 
asymmetry. The left hand side of Eq. (8) is Cy{pA), hence 
Eq. (8) can be rewritten as 
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Given the neutron-proton density profiles in the nucleus, 
from Eq. (9), a chosen value of 7 gives pA and hence 
CA- The one that satisfies Eq. (7) is the desired solution 
for 7. Once 7 is known, the equivalent density pA, the 
symmetry energy slope parameter and symmetry incom- 
pressibility are determined. From Eq. (7), there are two 
solutions for 7 and hence for pA- Calculations show that 
one of the solutions is unphysical as this gives pA ~ 5po- 
The procedure so described is expected to work best for 
heavy nuclei where volume effects predominate over those 
coming from surface. The heavy spherical nucleus ^°®Pb 
usually serves as a benchmark for extracting nuclear bulk 
properties, we choose this nucleus for our calculation. 

Apriori knowledge of Cy{pQ)^ Cs, po and the proton and 
neutron density distributions in the nucleus is required to 
extract values of 7 and p^- In different parametrizations 
of the nuclear masses |3ll - [35| . C„(po) is ~ 31 MeV, Cs 
is ~ 55 MeV; and po hovers around the canonical value 
~0.16 fm^'^ in different nuclear EOS models. As stated 
earlier, we take the values C^(po)=32.1 ±0.31 MeV and 
Cs =58.91 ±1.08 MeV. This value of C„(po) matches 
very well with 32.51 MeV, the one obtained from the 
latest mass systcmatics by MoUer et a/[l|; from their 
quoted value of 28.54 MeV for the surface stiffness pa- 
rameter Q, the value of Cs from ^"^Pb also comes very 
close, 58.16 MeV. For saturation density, we fix po =0.155 
± 0.008 fm^'^; this encompasses the saturation densities 
that come out from the EOS of different Skyrme and 
RMF models. The point proton distribution is known 
from experiments, the neutron density distribution is 
laced with much uncertainty though. 
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From a recent covariance analysis [25| , a lack of corre- 
lation of the neutron skin with some of the fundamental 
properties of nuclei like the isoscalar incompressibility, 
saturation density and the nucleon effective mass is sug- 
gested. The binding energy is also seen to be a poor 
isovector indicator. This is in consonance with the sug- 
gestion of effective nucleon-nucleon interactions of dif- 



(RMF), that give in the framework of niicroseopie mean- 
field theory different values of the neutron skin in ^°^Pb 
[21I l36j without compromising the basic nuclear prop- 
erties mentioned earlier. Aided by the further informa- 
tion that the neutron skin calculated with an effective 
interaction is strongly correlated with the corresponding 
symmetry energy slope parameter L , with empirical 
knowledge of Cv, Cs, po and the proton density distribu- 
tion, we now show how using Eqs. (7) and (9) for a heavy 
nucleus, both L and Rskin can be calculated. 

The parameters of the interactions BSR8-BSR14 (36| . 
FSUGOLD HI, NL3 ^ and TMl H] have been used 
to generate the proton and neutron density profiles of 
208pj^ in the RMF model. Alongwith many experimen- 
tally observed properties of finite nuclei and nuclear mat- 
ter, these interactions reproduce the proton r.m.s radius 
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Pb [Rp 
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=5.451 fm) extremely well. 



The neutron r.m.s radii vary considerably though, the 
calculated neutron skin varies from 0.17 fm to 0.28 fm. 
The symmetry energy slope parameter evaluated with 
these interactions using Eq. (1) are displayed in Fig. 1 
as a function of the corresponding calculated neutron 
skin Rskin (blue filled triangles). The magnitude of L 
increases with Rskin, its functional dependence L{Rskin) 
shows the usual linear correlation. These different in- 
teractions yield, using Eq. (1), different values of Cv{po) 
(31-38 MeV). We also use Eqs. (7) and (9) to calculate 
7 and hence L by employing the microscopic densities 
for the protons and neutrons obtained within the RMF 
models using the above mentioned parameter sets. These 
values of L are depicted in the same figure by the shaded 
region, the spread at a particular Rskin arising from the 
uncertainties in the values of C„, Cg and po. 

The shaded region projects a strikingly different de- 
pendence of L on Rskin- The filled red squares in the 
shade represent the median values for L, the filled green 
circles represent its lower as well as the upper bounds. 
The slope parameter L is seen to decrease here weakly 
with the neutron skin. This possibly originates from the 
fact that with given values of po, C„(po) and Cs, for a 
particular nucleus, Cv{po){^)'^ = asym{A) (Eq. (2)) is a 
fixed quantity, (pa)'^ is thus a constant for all the chosen 
interactions. Because pA and 7 have to satisfy Eq. (7), 
there is not much latitude in their values, resulting in the 
weak variation in L as shown. The intersection of the 
linear function L(Rskin) with the shaded region projects 
out those values of neutron skin for 208p|^ ^^lat are com- 
mensurate with the given "experimental" windows for 
Cv{po)tCs and pq. The corresponding calculated values 
of L are also accordingly projected out. The section de- 
picting the bounds of L and also of the neuron skin of 
208pj^ is shown by the box (magenta) in the figure. We 
find them to be L = 64 ±5 McV and i?sfc,;„(^°^Pb) = 
0.188 ± 0.014 fm. Our scheme for finding L is found to 
be quite robust. This is tested by choosing Woods-Saxon 
density profiles which are realistic but may not be very 
accurate. The proton density is adjusted to reproduce 
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FIG. 1: (Color online) The blue triangles represent L calcu- 
lated using Eq. (1) with different RMF interactions. They 
are plotted as a function of corresponding Rskin for ^"^Pb. 
The shaded region represents the envelope of possible L- values 
with different RMF interactions obtained using Eqs. (7) and 
(9). The acceptable window for the values of L and Rskin for 
^''*Pb is represented by the magenta box. 



tron density profiles are varied so that the entire range 
of 0.13 fm to 0.47 fm for Rskin as obtained from the 
PREX experiment is covered. Even in this large range 
of Rskin, L is confined within 55 to 85 MeV; the median 
value (70 MeV) is not much different from the one ob- 
tained with microscopic densities. Knowledge of L helps 
in predicting neutron skins of nuclei. Estimates for the 
neutron skin for a few nuclei are displayed in Fig. 2; they 
are obtained from the intersection of the shaded region 
showing the calculated limits of L =64 ±5 obtained from 
microscopic mean-field densities with the linear function 
L{Rskin) calculated for those nuclei in the RMF model 
with different energy density functionals. The values of 
the neutron skins displayed for the nuclei ^^^Sn, ^°Zr 
and ^^Ca are seen to be 0.196±0.014, 0.107±0.007 and 
0.182±0.008 fm, respectively. 

In summary, based on the empirical knowledge of the 
volume and surface symmetry coefficients and the nu- 
clear saturation density, we have presented a model that 
yields the density dependence of the symmetry energy in 
tighter bounds. This helps in making a precision pre- 
diction of the neutron skin of different nuclei, including 
the currently experimentally studied nucleus ^°^Pb, in 
PREX. We suggest that our determination on the lim- 
its on L and Rskin of experimentally studied nuclei be 
considered in properly evaluating nuclear energy density 
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FIG. 2: (Color online) The values of L are shown as a function 
of Rskin for ^^"^Sn (blue stars), ^'^Zr (magenta triangles) and 
''^Ca (green squares) calculated with different RMF interac- 
tions. The shaded region corresponds to the projected-out 
range in the values of L as obtained from Fig. 1. 
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